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Abstract. Delayed feedback laser dynamics is described by means of Lang-Kobayashi equation model. 
Since a lot of initial states asymptotically approach to periodic attractor in the phase space, only periodic 
steady-state regimes have been studied here. Lyapunov transformation allows us to reduce problem to 
the differential equation of the first order whereas the spectrum of laser oscillation is governed by the 
appropriate eigen value problem. Using the symmetry proves that there is the dual laser system with 
anticipated feedback which has the same dynamic characteristics as laser with delayed feedback. 
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1 Introduction 

A delayed optical feedback (DOFB) semiconductor laser is 
characterized by a variety of dynamical regimes pL| : some 
of them being used in optical communication systems. The 
complexity of these regimes is caused by the presence of 
two coupled resonators, one of which is formed by the 
reflecting end facets of a semiconductor crystal, and the 
other - by a combination of a highly reflecting end facet of 
the crystal and an external mirror. Lang and Kobayashi [2\ 
derived dynamic equations for a laser diode with an addi- 
tional external mirror by neglecting the reflection of exter- 
nal radiation from the diode facet. The Lang - Kobayashi 
(LK) equations form the basis of modern theories describ- 
ing the dynamics of DOFB laser diodes. The system of 
LK equations contains six key parameters, changing of 
which results in origin of a variety of dynamical regimes 
up to chaos. There are two alternative methods to pro- 
duce a variety of dynamical regimes in the semiconductor 
laser system without external mirror. In the first of them 
an optically injected diode laser is used. In the sec- 
ond one a system of the mutual optically coupled diode 
lasers [5J[Bj,l7 is applied. There is a common element which 
has all of these lasers. The part of the radiation is injected 
into laser resonator from the outside. It may be harmonic 
radiation as in OS], or undefined field in the general case. 

From mentioned point of view a delayed feedback laser 
and a system of the mutual optically coupled lasers are al- 
most equivalent. The reflection of the field from the mirror 
can be interpreted as an existence of the some radiating el- 
ement behind the mirror. So some virtual laser placed in 
image space exchanges radiation with the original laser. 
The characteristics of virtual laser and mutual coupling 
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constants are unknown and must be defined with the help 
of the theory. This problem is relatively simple if the laser 
regime is stationary. Yet it is well known that between 
this state and chaos, the delayed feedback laser demon- 
strates a lot of dynamical regimes. All of them belong 
to either transient states or to steady-state periodic [9] 
(conditionally-periodic [10]) oscillations. The existence of 
nontrivial dynamical properties of the nonlinear system 
is caused by delayed feedback. The new dynamic state of 
the system is appeared at a bifurcation points, sequence 
of which is developed with feedback variation. Asymptotic 
behavior of the dynamic state of the laser system is one of 
the steady-state regular oscillations, if the regime is not 
chaotic. The steady-state oscillations are reached in the 
inversion relaxation time. That is why except the short 
transient states the solving of the nonlinear system can 
be studied as a series expansion of the periodic function. 
It is natural that a numerical integration of LK equations 
is simple enough in comparison with a solving of this ex- 
pansion. But when DOFB laser can have several solutions, 
an infinite set of initial conditions must be verified to find 
some of them. Conversely, the our method can be repre- 
sented as a problem of eigen values, each being arranged 
by the energy stored in oscillations. To solve the nonlin- 
ear LK equations along periodic path it is used Lyapunov's 
transform [TT] . For that it is necessary to linearize system 
of equations with the help of small variations procedure 
along the closed path. Separating the tangent variations 
we find a new nonlinear equation with the same temporal 
evolution as initial problem. By applying mirror symme- 
try a theory permits to find parameters of virtual laser 
located in the image space. The dual laser in image space 
has anticipated dynamic with respect to original. So used 
concept couples a dynamics of DOFB laser with proper 
dynamics of the mutual optically coupled lasers. 
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2 Model of nonlinear laser dynamcis 

The nonlinear dynamics of the semiconductor laser with 
delayed optical feedback (DOFB) arising from the reflec- 
tion of the part of the radiation from external mirror is 
described by Lang-Kobayashi (LK) equations [2]: 

— = (1 - iR)NE + iMe'^'Eit -td), 

(1) 

dN 9 
T— = P -N -(1 + 2iV)|£;r. 
dt 

The first equation is used for description of the slowly 
varying envelope of the field amplitude E{t), the second 
one describes the inversion dynamics N(t) in relation to 
dimensionless time. R is the linewidth enhancement factor 
being responsible for reflective index changes with increas- 
ing gain. The module and phase of the feedback is denoted 
as M and k accordingly. The time lag of the argument in 
the delayed term of the equation is defined by the transit 
time over external resonator td ■ The field gain is propor- 
tional to population inversion, the dynamics of which is 
described by the second LK equation [SJ|H] . Dimensionless 
value N is expressed by means of carrier concentration Nc 
over threshold Nth'- N = l/2gTph{Nc — Nth), where g = 
{c/n)dG/dNc is a differential gain of the active medium. 
Here c and n are the speed of light and group refraction in- 
dex. The photon lifetime is t~i^ — {c/n){L~^\nr~^ +cec), 
while Qfc is distributed losses in passive structures, diode 
cavity length equals L and amplitude reflectivity of the 
end facets corresponds to r. The threshold carrier density 
Nth = Ntr + (gTph)^^ takes into account active losses and 
passive losses of the radiation inside volume and passive 
losses on the chip facets. The carrier density for trans- 
parency Ntr is used for denotation of the active losses. 
Here P — l/2gTph{jTs — Nth) is normalized pump inten- 
sity, and Ts is a carrier lifetime without induced transi- 
tions, J is a constant rate at which carriers are injected 
into active region. Dimensionless time is implemented with 
the help of the normalization to value photon lifetime Tph, 
which accounts for both passive in-cavity losses and losses 
on diode facets. In particular, the dimensionless inversion 
relaxation time is defined as: T = Ts/Tph- The typical value 
for this parameter is TsalOOO. The intensity of the field can 
be evaluated from the dimensionless E as a, photon flux 
per unit area 2cE^ /ngTs with energy hu! per photon. 

The LK equations system is provided with memory 
about previous states. As a result wide variety of differ- 
ent states [T] can be found as a function of initial distri- 
bution of the field and population inversion over time in- 
terval equaling time delay. Nevertheless, when pump level 
is small enough, randomly chosen set of initial conditions 
lead to regular dynamics. Based on theoretical grounds of 
the Kolmogorov-Arnold-Moser theory [T^] if the perturba- 
tion formed by delayed feedback is small enough, most of 
the phase trajectories of the system belong to the periodic 
invariant torus. Therefore it is very important to inves- 
tigate the steady-state periodic oscillations for LK equa- 
tions. The mathematical aspects of the theory of equations 
with delayed feedback are under construction |13) . 



In this work the regular periodic steady-state solutions 
of the LK equations are studied with the help of some 
methods of the linear theory of perturbations. When pe- 
riodic oscillations are steady-state, the perturbations are 
described by linear differential equations with periodic co- 
efficient. The system is complicated, but the certain equiv- 
alent system ordinary equations can be found from the 
fact of solution periodicity. The reducibility transforma- 
tion was proved by Lyapunov by using the property of 
the solution periodicity only '11'. This method helps us 
to consider the delayed feedback for periodic steady-state 
modes as some set of the effective DOFB constants. The 
combination of linear problem for perturbations and equa- 
tions for effective coefficients of DOFB permits us to solve 
the inverse problem of finding nonlinear oscillations for LK 
equations. 

After substitution E ~ Eg exp{i/3t + ip{t)) the dynamic 
of the field can be studied with the help of complex phase 
function ip. Unknown parameters were eliminated in ac- 
cordance with stationary state: ^ = 0, ^ = i/3, where 
j3 = M(cosx + Rsinx) and x = P'''d — k which is a new 
representation of the feedback phase. The stationary value 
of the population inversion Nq = —Msinx defines the 
filed intensity E^ ^ {P - No)/{l + 2No) in the station- 
ary regime, and /3 is a frequency detuning from the center 
of the reference wave. If the feedback strength is small 
enough (and hence Nq <^ 1), the equations for variables 
and n = N — Nq are: 

= (1 - iR)n + iMe-'"" [exp(V'(t - td) - Tpit)) - 1] 

(2) 

Here we have used the following definitions: = 2{P — 
iVo)/T w 2P/r and Tr^ = T-\l+2P)/{l+2No) ^p-^. 
To study the behavior of the perturbations to the solution 
of the system ([2]) it is important to take into considera- 
tion that usually tUr ^ T~^, that is, the characteristic 
frequency of oscillations is much higher than inverse time 
of the inversion relaxation. So our description is bounded 
only to those frequencies, which obey this condition. 



3 Perturbation analyze of LK equations 

To analyze the features of the equations ^ one can per- 
form a transfer to linear equations for small perturbations. 
With the help of substitution ^ tfj + Sip, n ^ n + 5n 
from initial nonlinear equations it is easy to produce lin- 
ear system of equations for small variances from exact 
solution. To study system of equation for variances from 
periodic solution, we can consider all coefficients of such 
system as being formed from periodic functions. With the 
help of Lypunov transformation X ~ L(t)Y this system 
X = P{t)X can be reduced to the system of the equation 
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with constant coefficients. In tfie one-dimensional case the 
reduced system is [TT] : 

Y = iujY. 

This equation corresponds 

S = X/iX = u} + L/iL, 

for original variable X. Because of the Lyapunov trans- 
formation L{t) is periodic function by definition, the de- 
formed limit cycle in this space is governed by the nonlin- 
ear equation: 



<P = U} 



Sfc exp {ik<P), 



where Sfe is some constant, obtained as result of series ex- 
pansion procedure for the periodic function L/L. So the 
Lyapunov transformation defines new nonlinear equation, 
which has the same temporal evolution as the original sys- 
tem. The next calculations are directed to determine the 
specific form of the expansion constants. 

Without DOFB the system of equations ^ has sta- 
ble stationary state ip = 0, so any initial perturbation 
decreases with time. DOFB changes the behavior of the 
system and this causes appreciable structural reconfigu- 
ration of the phase-plane portrait. The different types of 
the solutions are divided among them by the bifurcation 
points. 

To solve the linear system equations for small pertur- 
bations produced from ^ we divide the system into two 
different parts; the homogeneous part A and inhomoge- 
neous part B with delayed feedback. The first one doesn't 
have the feedback term and can be easily solved, the sec- 
ond part is governed only by DOFB. The strong recon- 
figuration of the phase-plane portrait shows us that any 
attempt to find the solution as a finite expansion in series 
of powers of feedback coefficient M is incorrect. Further- 
more the homogeneous part of the solution is vanishing 
value as opposed to a specific solution. The specific solu- 
tion of this system can be produced from the homogeneous 
part of the solution with the help of constant variation 
method [14] . Suppose that Sip = dx-\-i 6y, dip* — Sx — i Sy 
and = 5^{t - td) - H{t). = - td) - 5i)*{t), 

Ait) = M exp(-ix + i'{t - Td) - ilAt))^ then: 



<5V(t) \ ; „ ; / ^A{t) 5i 

5'iIj*{t) \ ^ dt r(exp(/ dt'A)) -iA*{t)SC 
6n{T) I -I { V 

(3) 

here r(exp) is a chronological exponent [TS], an external 
integral in formula Q is an indefinite one. The matrix A 
can be written as: 

1-iR' 

l + iR 

^(w2/2)e2" -(w2/2)e2« _i/T 

where the new symbolic designation have been introduced: 
2u= {ip + ^f). 



As the diagonal term in the matrix A is proportional 
to —1/T, the homogeneous part of the solution relaxes to 
zero with time. The only specific solution of the inhomoge- 
neous system defines the dynamics at relatively long time 
intervals. That is why the behavior close to periodic at- 
tractor is defined by the integration interval in the vicinity 
of upper limit, when t t. 

The features of the solution are governed by the spec- 
trum of the matrix A. It should be taken into account 
that only two eigen values have time dependence: Aq = 
0, Ai_2 = — 1/2T± iwrexp(-u). So three eigen vectors of 
the matrix A can be expressed: 

A \ ^{l-^Ry 

Wo = -H , «i,2 = (l + iR) 

In three dimensional complex eigen vector space of the 
matrix A it is possible to find two dimensional real sub- 
space vq, vi +V2 in which the third line {dn component) is 
zero if oscillation frequency is greater than relaxation rate 
oJr ^ T~^. The choice like this results from the absence 
of feedback via population inversion, so the third line of 
the perturbation vector is filled with zeros. With the help 
of two-dimensional subspace one can exclude dependence 
on variable Sn and decrease the dimension of the equation 
system. The exterior integral in matrix equation ([3]) im- 
plies integration in small vicinity of the variable r. There- 
fore the inner integral located under the exponent is taken 
on infinitely small interval (t, t). Every integration step h 
on this interval doesn't change vector vq , but changes vec- 
tors vi,V2- Nevertheless their sum vi + V2 is again the sum 
of eigen vectors for the next step of integration with an 
accuracy o{h). This specific property corresponds to adi- 
abatic approximation. If we define the primitive function 
(p for exp(M) function as: 



dt'exp(u)= / dt' exp((V'-t-V'*)/2) = ^'(r)-<?(t) (4) 



and retain the same order of accuracy to chronological 
operator in equation ([3]), then: 

r(exp {] dt' A)){vi + V2) = 
t 

If periodic oscillations are steady-state, the system of 
linear differential equations for small perturbations will 
be analogous to differential equations with periodic coef- 
ficients. In this case the system can be reducible from the 
viewpoint of Lyapunov [llj . In order to implement Lya- 
punov's transformation to the system of ordinary equa- 
tions, the expression Mexp(— -I- ^(r — to) — tp{T)) can 
be expanded in the exp(ikuj(l>) function series, where num- 
bers k are whole numbers. This system of the functions 
is complete system, since exp(ifcw<?) = [exp(icj^)]'^ forms 
power series, but it is not orthogonal. 



A{t) = M exp(-ix + ^l^{t - td) - iAt)) 

OO 

k— — oo 



(5) 
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hk - are the effective feedback coefficients for harmonics 
with the number k. So the two equations system can be 
developed from the matrix equation ([3]) (it is clearly that 
Sip = 5<P{Retp) since In(^) = Re^', SA = A{t)5^ oc S<P): 



(H{t) 

1 -iR\ , ^ , / +1 

1 + iR J ^ ■ ' ■ . ,, ^ \-i 
q+ = -ojJ2 [Khk + h*_k)/'A exp(ifccj<P(t)), 
q- = Rq+ -ujY. [Hhk - h*_^.)/2i] exp(ifcw<?(t)) 



= JdtS<P X 
cos(a;r^(T) - ojr'l'it)) + q- 



(6) 



4 Nonlinear solution and mirror symmetry 

The system of equations ^ is correct for any small varia- 
tion towards nonlinear solution. If we take into considera- 
tion the variations Sip which are tangent to the nonlinear 
solution, 5ip o: ip, then we can convert the equations © to 
obtain an assumed solution of nonlinear equations The 
motion along actual path produces the canonical transfor- 
mations preserving the invariants. According to Noether's 
theorem 12J a invariance relative to the translation along 
vector tangent to the cyclic phase path leads to the system 
has the first integral. Consequently perturbation 5ip for es- 
tablished states should be searched proportional to ■(/'• The 
equations © formally have the form 6il) — f{'P)S<P(Reip). 
Hence there is some nonlinear master equation ip = F(<P), 
from which it follows given equation for small variation 
while SF/S<P = /(^)- Furthermore the differentiation of 
this nonlinear equation produces also the new expression 
'ip = fi<l>)d> = f{<I>)<P{Re'ilj). The real part of the tangent 
variations HeSijj can be expressed in terms of function 
which was introduced in Value of the ^ defines the 
phase of the nonlinear oscillations and its time dependence 
obeys the differential equation of second order: 



Retp = ^/<P = - 2^ ^3 



(few) - uj; 



{hk — h*_i^)exp{ikuj<l>) 



Here the k numerates motion on different average fre- 
quencies kLo<l>. The number of elements under a summation 
sign specifies actual dimension of the motion equations 
system. The more the number of components is used, the 
more detailed nonlinear dynamics is described. The order 
of the equation ^ can be decreased in the following way: 
by multiplying into the function (p and performing the in- 
tegration. As a result the equation of the first order was 
derived (sj. = Sk(hk,h*_fP) = s*_^): 



<^ - ^ Sfc exp(ifcw«?) = E, 



Sk 



1 {kio){hk - h*Lk) 

2 {kuf - ujI 



(8) 



In this equation the integration constant E appeared. The 
new condition to exclude this constant ca n be added. From 
the system of equations ([2]) it follows that exp(2ReV') — 1 — 
with accuracy of the small parameter M/P. The con- 
stant E fulfils the role of the energy level, and its meaning 



can be found from the expression ([S]) with the help of con- 
dition (P^ = 1. The converted condition is E<P = 1. The 
period of the oscillation repetition can be determined via 
the method of the separation of variables: 



E +} Skexp{ikO) ] =uj/f2 = E (9) 



k 



So, Lyapunov's transformation is fulfilled with the help 
of the L/L function expansion. In the limit of 1/T 0, 
this transformation results in integrability of the system. 
The constant uj is an eigen number, and pulse repetition 

frequency equals [2 = = ijj<1>. 

Thus, the complex amplitude of oscillation tp is evalu- 
ated in terms of the integral: 



^P{T) 



' j dt ^ dfc exp{ikuj(l>{t)) , 



-<2hk + (hk 



h-k) 



(10) 



(11) 



To define sequence of the constants hk in expansion 
of the A{t) function we introduce a basis function set 
Pt.k = exp(ifcw<?(i)). The conclusion is used here that all 
of them have the basic period of oscillation 27r/f?, while k 
is the harmonic number. To build the search algorithm of 
the constants hk we apply the method of Moore-Penrose 
pseudo inverse matrix Then a vector with compo- 
nents hk equals to: 

h = (P+P)-i(P+/l(t)), P+ = (P*)^ (12) 

Vector of the coefficients h represented by the for- 
mula (I12p is the best approximation in terms of least- 
squares if a finite number of basis functions are used. The 
matrix (P'^F)mk is the nonsingular Hermitian matrix and 
its elements {2t: / Q)fk-m are specified by Fourier integral 
over period: 

fk-,n^{n/2n) J dteMMk-m)'P{t)) = f;^^^. 

The value of the matrix element is defined by differ- 
ence kernel, so any diagonal contains single element in all 
places. It was found that their meaning can be calculated 
directly from equation ([5]). If all values of the {hk ~ ^1^) 
are real numbers, the phase function has the reflection 
symmetry: (I>{t) = —<P{—t). In general case it is possible to 
construct the correspondence ^'(i) = —<!>{— t). This dual 
function is again the solution of the equation ^ , and yet 



the numbers 



Ski~h- 



-h^.) = s^. have been mod- 



ified. The dual system solution is spectral-symmetric to 
original one. In accordance with (I5ll2p effective feedback 
coefficients hn satisfy the following equation: 

hn = Afe^*^^ / dr (i>{T) exp {-inuj<l>{T)) x 
27r 



J dt ^dfcexp(ifcw<l'(t))j , (13) 
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here the inverse matrix factor (P+P)^^^ — (a;/27r)s„_fc 
was properly taken into account via (P and with the help 
of the equation Leading diagonal of this matrix is 
filled with numbers sq = E, in the other cases numbers 
Sfe which were defined earlier in formula ^ under k ^ 
are used. It is obvious that elements {2tt / Q)fk-n of the 
matrix (P+P)^^ can be produced by solving linear sys- 
tem: Skfk-n = (■f^/w)(Son, here Smn is the unit matrix and 
E^^ = n/oj. And finally coefficients dk = dk{hk, h*_f,) 
used in last equation (fT5)l are introduced in ([TlT) . 

To transform into dual system we symmetric reflect 
the index sign in (IT^ and the direction of the time axis 
should be changed too. As result the sign of the time delay 
td is changed: 

n/n 

h^n = Me"'^^ / dr (^'(r) exp {-inuj<P' (t)) x 
2vr J 

-7r/f2 

/ -y- \ 

exp —t / dt '^d^kCxp{ikoj<P'{t)) (14) 

So the transformation —h^„ ^ h'n and x — > x' (together 
with ~d^k ~> d'k = dk{h'k,h'*_f,), s'k = Sk{h' k,h' *_,.)) 
sets up a correspondence with the vector of the dual sys- 
tem. Then dual system is governed by the same equa- 
tion (jl3p only with opposite sign of the delay time —td- 
Having compared the expressions (jl3l I14p . the necessary 
condition for the feedback coefficient phase x' ~ X = oi' 
K + k' = 2 arctan(l/i?) can be found. If this condition is 
realized, the dynamics of the both systems have coincident 
period of oscillation 2-k jQ. It follows from the mirror sym- 
metry which is defined here as the time axis reflection and 
mirror inversion of the spectral components — /i-„ — > 
So both systems have the same eigen number w, but feed- 
back delay time transforms to feedback anticipate time. 
Therefore the effect of the delayed feedback can be inter- 
preted as the mutual optical coupling between the laser 
and its virtual image. However coefficients of the mutual 
optical coupling have difference in phase. So, the dynamics 
of the laser with delayed optical feedback and response of 
proper pair of the optical coupled lasers are equiv- 

alent, while the regimes of the periodic oscillations are 
used. The time lag, equaled to time delay for DOFB laser, 
retains between coupled lasers with dual dynamics. 

d 

= (1 - iR)n + iMe^'^ [cxp(V'' - V) - 1] 

dt 

(15) 

d 

= (1 - iR)n - iMe^'^ [exp(V' - ip') - 1] 



5 Conclusions 

For LK equations we have developed new spectral method 
of the solving of the nonlinear dynamic states. Without 
direct integration of original LK equations it is possible 
to analyze the periodic steady-state modes of diode laser 



with delayed feedback. It was found that the periodic at- 
tractor evolution is described with the first order differ- 
ential equation. This description was introduced with the 
help of Lyapunov's theoretical concepts about reducibil- 
ity of the linear differential equations with periodic coeffi- 
cients to ordinary one. The constants in this equation are 
defined via spectrum decomposition of the effective de- 
layed feedback. Calculation of these coefficients is a non- 
linear algebraic problem. The algebraic equations are de- 
veloped via the Fourier transform with an anharmonic 
oscillation phase, which obeys the above-mentioned dif- 
ferential equation of the first order. Finally it was shown 
that for any DOFB laser the dual laser system with antic- 
ipated feedback can be found with the help of symmetry 
principle. The optically coupled laser and its dual image 
produce the same dynamics as a single DOFB laser. 
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